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Abstract. - In the framework of intermediate wave-packets for treating flavor oscillations, 
we quantify the modifications which appear when we assume a strictly peaked momentum 
^ ■ distribution and consider the second-order corrections in a power series expansion of the energy, 

(f) ' By following a sequence of analytic approximations, we point out that an extra time-dependent 

phase is merely the residue of second-order corrections. Such phase effects are usually ignored 
in the relativistic wave-packet treatment, but they do not vanish non-relativistically and can 
introduce some small modifications to the oscillation pattern even in the ultra-relativistic limit. 
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Over recent years, the quantum mechanics of oscillations [1-3] has experienced much 
progress on the theoretical front [4], in particular, not only in phenomcnological pursuit of a 
more refined flavor conversion formula [5-7] which, sometimes, deserves a special attention, 
• but also in efforts to give the theory a formal structure within quantum field formalism [8,9,11]. 
Under the point of view of a first quantized theory, the flavor oscillation phenomena discussed 
in terms of the intermediate wave-packet approach [12] eliminates the most controversial points 
rising up with the standard plane-wave formalism [13,14]. In fact, wave-packets describing 
propagating mass-eigenstates guarantees the existence of a coherence length [12], avoids the 
ambiguous approximations in the plane wave derivation of the phase difference [15] and, un- 
der particular conditions of minimal slippage, recovers the oscillation probability given by the 
standard plane wave treatment. Otherwise, strictly speaking, the intermediate wave-packet 
formalism can also be refuted, for example, in the context of neutrino oscillation since such 
oscillating particles are neither prepared nor observed [4] in this case. Some authors sug- 
gest the calculation of a transition probability between the observable particles involved in 
the production and detection process in the so-called external wave-packet approach [4,9,16]: 
the oscillating particle, described as an internal line of a Feynman diagram by a relativistic 
mixed scalar propagator, propagates between the source and target (external) particles rep- 
resented by wave packets. It can be demonstrated [4], however, that the overlap function of 
the incoming and outgoing wave-packets in the external wave-packet model is mathematically 
equivalent to the wave function of the propagating mass-eigenstate in the intermediate wave- 
packet formalism. Thus, as a preliminary investigation concerning with the existence of an 
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extra time-dependent phase added to the standard oscillation term a 2 ™ q * [14], we avoid the 
field theoretical methods in detriment to a clearer treatment with intermediate wave-packets 
which commonly simplifies the understanding of physical aspects going with the oscillation 
phenomena. 

The main aspects of oscillation phenomena can be understood by studying the two flavor 
problem. In addition, substantial mathematical simplifications result from the assumption 
that the space dependence of wave functions is one-dimensional (z-axis). Therefore, we shall 
use these simplifications to calculate the oscillation probabilities. In this context, the time 
evolution of flavor wave-packets can be described by 

$>(z,t) — <fii(z, t) cos 9 i> x + <p 2 (z, t) s\n9 v 2 

= \4>i{z, t) cos 2 9 + <j) 2 (z, t) sin 2 9] v a + [<t>i{z, t) — (fi 2 (z, t)] cos sin vp 

= <f> a (z,t;0)i/ a + <f>p(z,t;9)i>f3, (1) 

where u a and vp are flavor-eigenstates and u : and u 2 are mass-eigenstates. The probability 
of finding a flavor state up at the instant t is equal to the integrated squared modulus of the 
vp coefficient 



P{v a -> vp-.t) = I dz 

•'-co 

where Fo(£) represents the interference term given by 



sin 2 [26»] r , . . . 

= ^ L ^{l-Fo(i)}, 



Fo(t) = Re 



dz 4>\{z,t) (f) 2 {z,t) 



(2) 



(3) 



Let us consider mass-eigenstate wave-packets given by 



0«(*,O) 



2 
ira 2 



exp 



exp [ipi z] , 



(4) 



at time t = 0, where i = 1, 2. The wave functions which describe their time evolution are 



.(*,*) = / 



dy z 



-iE$t + ip z z 



(5) 
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. In order to obtain the 



where EpJ = (pi +m 2 i ) 2 and <p(p z -Pi) = (2ira 2 ) 4 exp 

oscillation probability we can calculate the interference term Fo(f) by solving the following 
integral 



/ + oo 
<p(Pz-Pi)<p(Pz-P2) exp [-i AE (Pz) t] = 



exp 



-(aAp) 2 



dp z 
2tt 



f 2 (Pz-Po) exp [-i AE (Pz) t], 



(6) 



where we have changed the z-integration into a p^-integration and introduced the quantities 
Ap = p!-p 2 , p„ = + p 2 ) and AE Pz = Ep] -Ep 2 } The oscillation term is bounded by the 
exponential function of a Ap at any instant of time. Under this condition we could never 
observe a pure flavor-eigenstate. Besides, oscillations are considerably suppressed if a Ap > 1. 
A necessary condition to observe oscillations is that aAp <C 1. This constraint can also be 
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expressed by Sp 3> Ap where Sp is the momentum uncertainty of the particle. The overlap 
between the momentum distributions is indeed relevant only for Sp 3> Ap. Consequently, 
without loss of generality, we can assume 



Fo(t) 




-p ) exp [-i AE Pz t] 



(7) 



In the literature, this equation is often obtained by assuming two mass eigenstate wave- 
packets described by the "same" momentum distribution centered around the average mo- 
mentum p — p . This simplifying hypothesis also guarantees instantaneous creation of a pure 
flavor eigenstate v a at t = [15]. In fact, for <j>i(z, 0) = (f> 2 (z,0) we get from Eq. (JIJ 



2 



exp 



exp [ip z] 



(8) 



and (j)p(z,Q,6) = 0. In order to obtain an expression for cf> i (z,t) by analytically solving the 
integral in Eq. (JjjJ we firstly rewrite the energy EpJ as 



E 



W _ 



E, 



Ef 



Et [l + <Ti (a, + 2v,)p 



(9) 



where E t = (m 2 + p 2 )? , v 4 = and <r 4 = p ^ Po . The use of free gaussian wave packets 
[4, 5, 9, 10] is justified in non-relativistic quantum mechanics because, in most of the cases, 
the calculations can be carried out exactly for these particular functions. The reason lies in 
the fact that the frequency components of the mass-eigenstate wave-packets, Ep} — pl/2m i , 
modify the momentum distribution into "generalized" gaussian, easily integrated by well 
known methods of analysis. The term p\ in Ep'J is then responsible for the variation in time 
of the width of the mass-eigenstate wave-packets, the so-called spreading phenomenon. In 
relativistic quantum mechanics the frequency components of the mass-eigenstate wave-packets, 



do not permit an immediate analytic integration. This difficulty, however, 



= VpT 

may be remedied by assuming a sharply peaked momentum distribution, i. e. (aE^ 1 ~ 
(Tj <C 1. Meanwhile, the integral in Eq. (JHJ can be analytically solved only if we consider terms 
up to order o~ in the series expansion. In this case, we can conveniently truncate the power 
series 



Pz 



E, 



Et+p^ + ^a 2 . 



(10) 



and get an analytic expression for the oscillation probability. The zero-order term in the 
previous expansion, E tl gives the standard plane wave oscillation phase. The first-order term, 
PoCTi, will be responsible for the slippage due to the different group velocities of the mass- 
eigenstate wave-packets and represents a linear correction to the standard oscillation phase 

2 

[15]. Finally, the second-order term, ^-cf, which is a (quadratic) secondary correction will 
give the well-known spreading effects in the time propagation of the wave-packet and will be 
also responsible for a new additional phase to be computed in the final calculation. In the case 
of gaussian momentum distributions for the mass-eigenstate wave-packets, these terms can 
all be analytically quantified. By substituting 1)10(1 in Eq. © and changing the ^-integration 
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into a di-integration, we obtain the explicit form of the mass-eigenstate wave-packet time 
evolution, 



where 8 t (t, z) = || 



7raf(i) 



arctan 



2 m t 



exp [— i (0j(f, z) + Eit — pa z)] exp 
| and fli(t) = a 



f) 2 



« 2 W 



(11) 
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4m; 



. The time- 



2mft (z-Vj t) 2 
d^Ef af(t) 

dependent quantities a ; (t) and 9i(t, z) contain all the physically significant information which 
arise from the second order term in the power series expansion 1)10(1. By solving the integral 
(|7| with the approximation and performing some mathematical manipulations, we obtain 

Fo(f) = Bnd(£) x Osc(t), (12) 

where we have factored the time- vanishing bound of the interference term given by 

(Avf) 2 



Bnd(£) = [1 + Sp 2 (f)] 3 exp 



2a 2 [1 + Sp 2 (f)]_ 

and the time-oscillating character of the flavor conversion formula given by 

Osc(f) = Re{exp[-iAEt - i@(t)}} 
= cos [A.Ei + 0(f)], 



where 



and 



Sp(t) = £A (£) 



= P 



AVt 



0(f) = \ arctan [Sp(i)] 



2 2 
" /'() 
~2p T ~ 



, Sp 3 ( t) 
[i+Sp 2 ( t )j 



(13) 

(14) 
(15) 
(16) 



with p=l- 3 + (^f ) 2 f£ and £ = \jE~TEl. The time-dependent quantities Sp(t) and 0(f) 
carry the second order corrections and, consequently, the spreading effect to the oscillation 
probability formula. If AE -C -E, the parameter p is limited by the interval [1,-2] and it 

2 

assumes the zero value when sa -|. Therefore, by considering increasing values of p , from 
non-relativistic (NR) to ultra-relativistic (UR) propagation regimes, and fixing , the time 

2 

derivatives of Sp(f) and 0(f) have their signals inverted when -g^ reaches the value |. The 
slippage between the mass-eigenstate wave-packets is quantified by the vanishing behavior of 
BND(f). In order to compare BND(f) with the correspondent function without the second 
order corrections (without spreading), 



BND ws (f) = exp 



(Avf) 2 
2a 2 



we substitute Sp(£) given by the expression <|14|l in Eq. (|13fl and we obtain the ratio 



BND(f) 
BND ws (f) 



AEt 

a 2 E 2 



exp 



p 2 p 2 (AEt) 4 



2 a 6 E s 



1 



' AEt \ 2 
\a 2 E 2 ) 



(17) 



(18) 



The NR limit is obtained by setting p 2 = 1 and p„ = in Eq. 11711 . In the same way, the UR 
limit is obtained by setting p 2 — 4 and p a — E. In fact, the minimal influence due to second 
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order corrections occurs when J=j^ ss i (p«0). Returning to the exponential term of Eq. (JTSJl, 
we observe that the oscillation amplitude is more relevant when Avt -C a. It characterizes the 
minimal slippage between the mass-eigenstate wave-packets which occur when the complete 
spatial intersection between themselves starts to diminish during the time evolution. Anyway, 

under minimal slippaqe conditions, we always 

have BNDW ; „ L 

DNDws(t) 

The oscillating function Osc(i) of the interference term Fo(t) differs from the standard 
oscillating term, cos [AEt], by the presence of the additional phase Q(t) which is essentially 
a second order correction. The modifications introduced by the additional phase &(t) are 
discussed in Fig. ^ where we have compared the time-behavior of Osc(i) to cos[AEt] for 
different propagation regimes. The bound effective value assumed by O(i) is determined by 
the vanishing behavior of BND(t). To illustrate this flavor oscillation behavior, we plot both 
the curves representing Bnd(£) and 8(i) in Fig. [5] We note the phase slowly changing in the 

NR regime. The modulus of the phase |8(t)| rapidly reaches its upper limit when > | and, 
after a certain time, it continues to evolve approximately linearly in time. But, effectively, the 
oscillations rapidly vanishes. By superposing the effects of Bnd(£) in Fig.[3and the oscillating 
character Osc(i) expressed in Fig.^ we immediately obtain the flavor oscillation probability 
which is explicitly given by 



sin 2 [201 



1-[1 + Sp 2 (T)] 4 exp 



(Avt) 2 



2a 2 [1 + $P 2 (t)\ 



cos [AEt + Q(t)] 



(19) 

and illustrated by Fig. Obviously, the larger is the value of a E, the smaller are the wave- 
packet effects. If it was sufficiently large to not consider the second order corrections of Eq. JHJ, 
we could rewrite the probability only with the leading terms (slippage effect), 



P(v„ 



sin 2 [201 



1 — exp 



(Avt) 2 
2 a 2 



cos [AEt] 



(20) 



which corresponds to the same result obtained by [15]. By assuming an UR propagation 
regime with t ~ L and Ei ~ p , under minimal slippage conditions (Av L <C a), the Eq.J^OJ) 
reproduces the standard plane wave result, 



P(u a -> v$\ L) 



sin 2 [20] 





Am 2 




1 1 — COS 




} 




. 2p 





(21) 



since we have assumed a E ^> 1 . 

By summarizing, we have obtained an explicit expression for the flavor conversion formula 
for (U)R and NR propagation regimes which is valid under the particular assumption of a 
sharply peaked momentum distribution. We have also observed that the spreading represents 
a minor modification effect which is practically irrelevant for (ultra)relativistic propagating 
particles. In particular, the intermediate wave-packet prescription elaborated here can be 
discussed in the context of neutrino flavor oscillations. We have concentrated our arguments 
on the existence of an additional time-dependent phase in the oscillating term of the flavor 
conversion formula. Such an additional phase presents an analytic dependence on time which 
changes the oscillating character in a peculiar way. These modifications are minimal when 
pi w ^E 2 and more relevant for NR propagation regimes. The existence of an additional 
time-dependent phase in the oscillating term of the flavor conversion formula coupled with 
the modified spreading effect can represent some minor but accurate modifications to the 
(ultra)relativistic oscillation probability formula which leads to important corrections to the 
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phcnomcnological analysis for obtaining accurate ranges and limits for the neutrino oscillation 
parameters. The relevance of such second-order corrections depends essentially on the value 
of the product between the wave=packet width a and the averaged energy flux E which 
parameterize the power series expansion here proposed and quantified. 

Finally, we know the necessity of a more sophisticated approach is understood. It involves a 
field-theoretical treatment. Derivations of the oscillation formula resorting to field-theoretical 
methods are not very popular. They are thought to be very complicated and the existing 
quantum field computations of the oscillation formula do not agree in all respects [4] . The 
Blasone and Vitiello (BV) model [8, 11] to neutrino/particle mixing and oscillations seems 
to be the most distinguished trying to this aim. They have attempt to define a Fock space 
of weak eigenstates to derive a nonperturbative oscillation formula. Also with Dirac wave- 
packets, the flavor conversion formula can be reproduced [17] with the same mathematical 
structure as those obtained in the BV model [8, 11]. In fact, both frameworks deserve a 
more careful investigation since the numerous conceptual difficulties hidden in the quantum 
oscillation phenomena still represent an intriguing challenge for physicists. 

* * * 
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Fig. 1 - The time-behavior of Osc(t) compared with the standard plane-wave oscillation given by 
cos [Ai5t] for different propagation regimes. The additional phase Q(t) changes the oscillating char- 

acter after some time of propagation. The minimal deviation occurs for -g^ ~ § which is represented 
by a solid line superposing the plane-wave case. 
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Fig. 2 - The values assumed by Q(t) are effective while the interference term does not vanish. In 
the upper box we can observe the behavior of Bnd(£) which determines the limit values effectively 

„2 

assumed by O(t) for each propagation regime. For relativistic regimes with > ±, the function 
6(t) rapidly reaches its lower limit as we can observe in the small box above. We have used aE = 10. 
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Fig. 3 - Flavor conversion probability time dependence obtained with the introduction of second- 
order corrections in the series expansion of the energy for a strictly peaked momentum distribution 
(0(of )). By comparing with the PW predictions, depending on the propagation regime, the additional 
time- dependent phase A4>(t) = AEt + 0(t) produces a delay/advance in the local maxima of flavor 
detection. Phenomenologically, it can introduce small quantifiable deviations to the averaged detected 
values of neutrino oscillation parameters. Essentially, it depends on the product of the wave-packet 
width a by the averaged energy E. 



